Theory of Excitons in Insulating Cu-Oxide Plane 
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We use a local model to study the formation and the structure of the low energy charge transfer 
excitations in the insulating CU-O2 plane. The elementary excitation is a bound exciton of spin 
singlet, consisting of a Cu + and a neighboring spin singlet of Cu-0 holes. The exciton can move 
through the lattice freely without disturbing the antiferromagnetic spin background, in contrast to 
the single hole motion. There are four eigen-modes of excitons with different symmetry. The p- 
wave-like exciton has a large dispersion width. The s-wave-like exciton mixes with p-state at finite 
momentum, and its dipole transition intensity is strongly anisotropic. The model is in excellent 
agreement with the electron energy loss spectra in the insulating Sr2Cu02Cl2- 
PACS numbers: 71.35.Cc, 74.72.Jt 



I. INTRODUCTION 

There have been extensive efforts in recent several 
years on the electronic structure in the layered Cu-oxides, 
in the hope to reveal the mechanism for the high T c super- 
conductivity. It is now well established that the low en- 
ergy physics of the undoped cuprates is well described by 
the spin-1/2 Heisenberg model in a square lattice, lead- 
ing to an antiferromagnetic insulator. Two-dimensional 
t — J model has been proposed and argued to describe 
the low energy physics for the cuprates |l],|2| , and it seems 
the model explains many unusual properties of the high 
T c materials. 

In this paper we present a theory for the lower en- 
ergy charge transfer excitations in the insulating Cu- 
oxide plane. Our work was motivated by the angle- 
resolved electron energy loss spectra (EELS) in the insu- 
lating Sr2Cu02Cl2, where the optical allowed transition 
shows a large energy dispersion width 1.5 eV, and the 
optical forbidden transition was observed to be strongly 
anisotropic j|. 

We use a local model to analyze the formation and 
the structure of the charge transfer excitation between 
Cu-3d x 2_ y 2 and 0-2pa in the Cu-oxide planes. The el- 
ementary charge transfer excitation is a bound exciton 
of spin singlet, consisting of a Cu + (quasi-particle) and 
a neighboring spin singlet of Cu-0 holes. The exciton 
can move through the lattice freely without disturbing 
the antiferromagnetic (AF) spin background, in contrast 
to the single hole motion. There are four exciton modes 
related to the four-fold rotation symmetry in the lattice, 
whose hopping matrices determine the dynamics and the 
symmetry of the charge transfer excitations. The p-wave 
like exciton has a large energy dispersion width. The 
s-wave-like exciton mixes with the p-state at finite k, 
and its optical transition intensity is strongly anisotropic. 
The model is in excellent agreement with the recent elec- 
tron energy loss spectra in the insulating Sr2Cu02Cl2. 
The theory also predicts a lower energy optical forbid- 
den mode, which should be observable in luminescence 
experiment Q . Since the mass of the particle-hole pair 
is heavier than the particle or the hole in a usual semicon- 



ductor or in other band- like insulator materials ||, the 
anomalous behavior of the dispersion widths between the 
single hole and the exciton in cuprates provides further 
evidence for the strong correlation in these compounds. 
A part of the results in this paper was briefly reported 
in a joined article with the EELS experiments j^]. The 
paper is organized as follows. In section [n|, we exam- 
ine the formation and the structure of small excitons in 
the insulating cuprates. The dynamics of the excitons is 
derived in section III. The exciton eigen-modes and the 



relevance to the charge transfer gap and to the EELS are 
described in section IV , Several remarks and discussions 
are made in section VJ 



II. STRUCTURE OF SMALL EXCITONS 

We consider a Cu-oxide plane, where the vacuum is the 
filled shell states of Cu + and 2 ~. In the undoped case, 
there is one hole in average in each unit cell of Cu02 , and 
the hole primarily resides on the Cu site and is of d x 2_ y 2 
symmetry. Because of the strong on-site Coulomb repul- 
sion between the Cu holes, the ground state is a spin- 
1/2 antifcrromagnet of Cu 2+ . To study the low energy 
charge excitations of a particle-hole pair, we shall focus 
on the transition from the Cu-3d a .a_ 1/ 2 to 0-2pa states. 
This transition is of particular interest due to the rela- 
tively small atomic energy difference and due to the large 
orbital overlap between the two states hence the strong 
intensity in experiments. The measured optical gap j(| 
is related to this transition. 

Let us start with a local pair of particle and hole and 
examine the formation and the structure of the exciton. 
As illustrated in Fig. [I], the quasi-particle is at a vacant 
Cu-site (Cu + ), and the quasi-hole is at the O-site. There 
are two important physical effects for the pair. The first 
is the Coulomb attraction between the particle and the 
hole, which should be considerably large in the insulat- 
ing compounds where the screening effect is weak. It 
is this attraction to bind the pair to form an exciton. 
For simplicity we shall assume the attractive potential 
to be short-ranged, which is — E c < when the particle 
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and hole are nearest neighbors (n.n.), and is zero oth- 
erwise. We expect that the more realistic potentials will 
not change the basic physics as long as the n.n. attraction 
dominates. The second is the strong Cu-0 hybridization 
which binds a hole on each square of O-atoms to the cen- 
tral Cu-hole to form a spin singlet j| . Let us denote this 
binding energy to be E s . As a compromise of the above 
two effects, we expect a spin singlet with more weight 
of the O-hole on the site near the Cu + . We can make 
the above argument more quantitatively using variation 
approach. Let the O-hole state be 

* - aP-i/2 + P + Py/2 - P-y/2) , U) 

where P is the atomic O-hole orbital, whose sub-index 
denotes the atomic position relative to their central Cu- 
hole (lattice constant=l), and P_£/2 is the one close 
to the Cu + . The total binding energy of the exciton, 
Ef,, takes maximum when a 2 = (1 + l/v 7 ! + 3c)/2, with 
c = E 2 /(2E C - E s ) 2 . Typically, a 2 = 0.5 or 0.75 for 
E c / E s =0.5 or 1. From the above analysis we see that 
the exciton in cuprates is a complex object, consisting of 
a Cu + and a nearby Cu-0 singlet. The spin triplet exci- 
ton energy has a binding energy of E C1 smaller than the 
binding energy Ef, of the spin singlet exciton, and will 
not be considered in this paper. 

The parameters E c and E s in cuprates should be 
within certain regions so that the system is stable against 
the excitonic excitations. Let e p be the atomic energy dif- 
ference between the Cu- and O-hole states, the stability 
condition in the local limit is Ef, < e p . If we further con- 
sider an extra hole (primarily residing on O-site) in the 
half- filled insulating Cu-oxides, the model system should 
be stable against excitonic excitation (one of the Cu-hole 
moves to the oxygen). This implies 2Et < e p + E s . Our 
discussion in this paper will be limited to the above phys- 
ical parameter region. 

There are four spin singlet exciton states for a fixed 
quasi-particle at R, denoted by t(R), where r is the 
coordinator of the Cu-site in the singlet relative to R. 
r = x, — x, y, and — y. These will be briefly denoted by 
x, x, y, and y respectively (see Fig. |l]). 

III. DYNAMICS OF EXCITONS 

The exciton can move through the lattice due to the 
n.n. Cu-0 hopping t p d and the 0-0 direct hopping t pp . 
The most important feature of the exciton dynamics is 
that the exciton motion does not disturb the AF spin 
background. This is because each exciton involves two 
neighboring Cu-sites, both are spinless. As far as spins 
are concerned, the spin singlet exciton is similar to a pair 
of bound n.n. holes, whose motion does not disturb the 
AF spin background. This property should be compared 
with the single hole dynamics in Cu-oxide, where the AF 
spin background is strongly disturbed, and the hole dis- 
persion is given by the Cu-Cu spin exchange interaction 



J instead of the hole hopping integral t, as shown theo- 
retically for the t — J model and observed in the recent 
angle-resolved photoemission experiments B 

We now study the exciton motion quantitatively. A 
general hopping process of the exciton in the coordina- 
tor space may be represented by r(R) — > r'(R') with the 
hopping integral t TT >(R— R'). Because the AF spin back- 
ground is unchanged Js) , the exciton motion is equivalent 
to that for a free particle. Let \GS) be the AF ground 
state, and 7^(R) is an operator to create an exciton r(R) 
in the ground state, so that 7^ (R)|G5) is the state of ex- 
citon t(R) in an otherwise AF spin background. The 
exciton hopping is described by an effective Hamiltonian 
H e ff, which acts on the Hilbert space of the single exci- 
ton states, 

H efn t(R)\GS) = t TT ,(R--R>)vl,(R>)\GS). (2) 

R',T' 

It follows that 

%^ T t(k)T(k) 7 (k), (3) 

k 

with 7 t(k) = (71,71,7^, 7^). T(k) is a 4 by 4 matrix, 
which determines the exciton dynamics, and T rr <(k) = 
J2 r t TT > (r)e lk ' r /N, with N the number of Cu-sites. Note 
that k is in the whole first Brillouin zone although there 
are two sub-lattices in the ground state. The matrix is 
hermitian, T TT /(k) = T*, T (k). Because t TT '(r) is real, we 
have T TT /(k) = T T * T ,(-k). The Cu-0 plane has certain 
symmetries. Taking into account of the Cu and O-hole 
orbital symmetries, we have T rr <(k) = T ff i{— k) from 
the parity invariance, T xy -(k x ,k y ) = —T xy (k x ,—k y ) from 
the inversion symmetry with respect to the x-axis, and 
T xx (k x , k y ) — Tyy(~k y , k x ) due to the four-fold rotational 
invariance. 

To pursue the theory further, we estimate the hopping 
integral t TT i (r) within the atomic limit and treat t p d and 
t pp perturbatively. We shall consider the limit a = 1 in 
Eq.(l) for simplicity and for the purpose of illustration. 
Up to the third order in perturbation in t pd or t pp , there 
are only four non-zero independent integrals: 

^2 = t x y{x) = tppt pd /(2epE c ), 

h = -t xx (x) = tl d [e~ 1 - (e p + Upp)' 1 }, 

t i = t xx (0)=t 2 pd /(2e p ). (4) 

A diagrammatic illustration for these hopping processes 
are shown in Fig. 0. In Eq. ([|) , XJ PP is the on-site O-hole 
Coulomb repulsion. We have kept the third order term, 
t2, because E c is small compared to e p . From Eq. (|J), 
<3 > <4. Since t pp > for the O-oxygen hole hopping 
, t2, £3, and ti are all positive. Using these t's and the 
symmetries discussed above, we obtain 
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TOO 



/ a(k x ) b(k) 6(k) \ 

a*(jfe x ) _0 6*(k) 6*(k) 

6*(k) 6(k) a(fc y ) 

^ 6*(k) b(k) a*(k y ) 



(5) 



where b(k x ,k y ) = —b(k x ,—k v ), and 
a(k x ) = U-t 3 e lk *, 



b(k) =ti+t 2 (e 1 *- +e-*"). 



(6) 



IV. EXCITON EIGEN-MODES 

iJ e // can be diagonalized to obtain the exciton solu- 
tions. For k along [11] direction, the analytical solu- 
tions are particularly simple and are listed in Table I. At 
k=0, the four exciton eigenstates have well defined local 
symmetry (9| , describing the uniformly moving molecular 
states. They are illustrated in Fig. [|. The s-wave and 
d-wave states are given by 



I*) = (" 
\d) = (- 



II 11 

2' 2'~2' 2 

III 1 

2' 2'2'~2 



(7) 



where the four components represent the amplitudes of 
the excitons 7 X , 7s , 7 a , 7j/ respectively. The energies for 
s- and d-states are 2bo — ao, and — 2&o — ao respectively, 
where a Q and 6 are their values at k = 0, and a = — 
ta > 0, bo = t\+t2- The two p- wave states are degenerate 
and the energy is — ao. In the parameter region — &o < 
a < b and b > (suitable for Cu-oxide, see Fi g. H 
caption), the s-wave state has the highest energy [ |10[ 7 
The optical experiment measures the dipole active mode 
at fc=0. Because of the d-symmetry of the Cu-hole, both 
the s- and d-wave modes are optically forbidden. The p- 
wave modes are optically active, and the charge transfer 
gap is given by 



A 



E b - ao. 



(8) 



For finite k along an arbitrary direction k, the eigen- 
state of the exciton is a linear combination of the molec- 
ular states. We shall denote the four exciton modes by 
Pi, P%, S, and D, according to their symmetries p\ , pi, s, 
d x 2_ y 2 in the limit k — > 0. For the twop-states, we define 
a dipole active mode \pi) and a dipole inactive mode \p2) 
as 



\pi) = (-k x ,-k x ,k y ,k y )/(V2k), 
\p 2 ) = (k y ,k y ,k x ,k x )/(V2k). 



(9) 



The exciton dispersions are plotted in Fig. ^. The pa- 
rameters are chosen |ll| to fit the measured EELS ||. 
The dispersions are anisotropic. The E — k relation of 
the P\ mode is monotonic, and the energy is peaked at 



k = (tt, 7t). The dispersion width is 2t 3 ~ 2t^ d /e p , which 
is a large energy scale, order of the hole hopping integral 
in cuprates, but much larger than the dispersion width 
of the single hole motion in the AF background which 
is ~ J f|. This result is consistent with the EELS [||. 
The dispersion along [10] direction is much flatter, also 
consistent with the EELS. As we can see from Fig. || the 
lowest energy mode is the -D-mode. This is in agreement 
with the week coupling method of Littlewood, Varma, 
and Schmitt-Rink @. 

In the EELS, an incident electron is inelastically scat- 
tered to create a pair of electron-hole with energy and 
crystal momentum (u>, k). The EELS directly mea- 
sures Im[— l/e(k, lu)), with e(k, uS) the dielectric func- 
tion, hence probes the dispersion and the symmetry of 
the charge transfer excitation fl3| . For an exciton, the 
intensity in the imaginary part of e(k, u>) around a pole 
is, 



I = const * k 



r \GS)\ 



(10) 



where ^ ex is the exciton state. For small k, the most 
dipole active excitations in the EELS can be identified 
to be the P\ mode. The dispersion obtained from the 
EELS |l| compares well with the theory. As k increases, 
the S- and D-mode excitons gradually become EELS ac- 
tive because of the mixing of the molecular state and 
because of the quadrapole contribution. It is difficult to 
make a quantitative comparison between the dipole and 
quadrupole contributions in the S- or D-mode. However, 
if the radius of the Cu-3d state is much smaller than the 
lattice constant, the dipole contribution is more impor- 
tant. We expect this to be the case in cuprates. Within 
the dipole approximation, we have 



J = Jo|<*exbl>r 



ID 



where Iq is the dipole active intensity of Pi-mode at k=0. 
Note the total dipole active intensity of the four modes 
is conserved to be Iq. In Fig. |[ we show the calculated 
intensity I of the excitons in Eq. (|ll]) as functions of 
k. Along [10] direction, as k increases, the intensity of 
the S'-mode increases rapidly, while the intensity of the 
Pi-mode decreases. The intensity along [11] is quite flat 
for the Pi-mode and remains to be zero for the 5-mode. 
This anisotropy is due to the symmetry of the exciton 
states. The S*-mode eigenstate along [11] contains only 
s and P2 states, both are dipole inactive. Furthermore, 
the transition matrix from the quadrupole contribution 
is proportional to (^ea; |e lk ' r |G5') , which is (k x — ky) for 
the s-state [[l4| . Therefore, the intensity for the >S-mode 
along [11] vanishes up to the quadrupole order. These 
features are in excellent agreement with the EELS ex- 
periment ||, where an optical forbidden transition was 
observed along [10] but not along [11] direction, and the 
intensity of the optical active transition along [10] di- 
rection is observed in the experiment to decrease as k 
increase. 
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The D-mode exciton, expected from the theory, is an- 
other optical forbidden state. The current EELS has not 
revealed this mode, probably due to the limit of the ex- 
periment resolution. Further spectroscopy measurements 
are needed to verify this rf-wave-like state. This mode 
can be active in the phonon-assisted optical process, and 
should be observable in the luminescence experiment. A 
recent optical absorption measurement indicates a very 
weak absorption state at about 0.5 eV in the undoped 
cuprates [|l4| . The d-symmetry state is an alternative to 
the magnon state or crystal field exciton proposed ear- 
lier for this weak absorption where the phonons could be 
involved. 



for the insulating Cu-oxides. The exciton moves through 
the lattice almost freely in the AF spin background. The 
model is in excellent agreement with the recent EELS in 
the insulating Si^CuC^C^. 
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V. DISCUSSIONS 



In this paper we discussed small excitons in the in- 
sulating Cu-oxide planes. These are the charge transfer 
excitations, which are separated from the lower energy 
excitations by an energy order of e p . The lower energy ex- 
citations are the magnons described by spin -1/2 Heisen- 
berg model at half-filled. In the doped case, the lower 
energy charge excitations are given by the hole motion in 
the AF spin background described by the 2-dimensional 
t — J model. 

We have used a local model to describe the excitons 
in Cu-oxide planes. If we include the individual kinetic 
energies of the quasi-particle and the quasi-hole, the ex- 
citons should have spatial extensions. The size of the 
exciton is determined by the balance of the kinetic ener- 
gies and the Coulomb attraction E c . This spatial exten- 
sion will lower the exciton energy. The exciton motion 
will be more complicated, and will contain an incoherent 
part disturbing the AF background. The large disper- 
sion width observed in the EELS may be viewed as an 
experimental indication that the spatial extension of the 
exciton does not play an important role and the small 
exciton proposed here may be a good approximation. 

We have so far not considered the electron-hole (e-h) 
continuum, which starts at energy A e _^ = e p — E s + 
Exin , where Exin < is the lowest kinetic energy of the 
independent Cu + and the formal Cu 3+ . For the cuprates, 
I-Eifini ^1 — 2 eV. We expect the e-h continuum to start 
above the k=0 p\ mode, and the exciton spectra extend 
into the e-h continuum due to the large spectral disper- 
sion. The states within the continuum region will then 
be damped, but can exist as resonant states contributing 
to the EELS. Since the contributions to the EELS from 
the e-h continuum have much less k-dependence, the ex- 
citons become the dominant source of the k-dependent 
spectra in the EELS. 

Our model applies to the insulating phase. In the 
metallic phase the Coulomb attraction between the quasi- 
hole and quasi-particle is substantially reduced due to the 
metallic screening, and is too weak to bind a pair. 

In conclusion, we have studied a local exciton model 
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FIG. 1. Structure of local excitons in insulating Cu-O 
planes. The open circles represent O-atoms, the solid circles 
represent Cu-atoms. (a) exciton x(R): the quasi-particle is at 
the vacant Cu-site R, and the quasi-hole is on the square of 
O-atoms with more weight at the left atom and forms a spin 
singlet with the central Cu-hole at position R + f. Not shown 
are a single hole on all other Cu-sites. (b) a simpler repre- 
sentation of (a), (c), excitons ic(R), J/(R) and y(R), with the 
spin singlet at different positions relative to R. 

FIG. 2. A list of five most important hopping processes 
from x(R) to x(R') and y(Rf). The other hoppings can be 
obtained by symmetry. The corresponding effective hopping 
integrals are given by the Eq. (0). 



Mode 


energy eigenstatc 


symmetry at k — > 


Pi 


-b+\a-b\ (-r?-, -1/2,7?-, 1/2) 


Pi 


P 2 


b-\a + b\ (-77+, 1/2, -77+, 1/2) 


P2 


S 


fe+|a + 6| (77+, 1/2, 77+, 1/2) 


s 


D 


-b-\a-b\ (ij_,-1/2,-»j_,1/2) 





TABLE I. Solutions of the four exciton modes for k along 
[11] direction. The k-dependence in a, b and b is implied, see 
Eq.(g|). I7± = |(a±6)/|a±6|. The last column applies to the 
region (—60 < «o < &o) suitable to cuprates, and ao, 60 are 
their values at k — 0. 



FIG. 3. Symmetries of the four molecular exciton wave- 
functions, s- and d- waves in Eq. (Q), and the two p- waves in 
Eq. (|^) with kx = k y . The center d- wave represents Cu + , 
a vacancy at Cu-site. The oxygen 2p x (2p y ) hole wavefunc- 
tions are represented by their orbitals. Only charge degrees 
of freedom are shown. The implicit spins are the same as in 
Fig. § 



FIG. 4. Calculated exciton energies E (eV) as a function 
of k= (k x ,k y ), for the parameters t\ = 0.4, ti = 0.126, 
£3 = 0.85, <4 = 0.65 eV. The energy is measured relative 
to the pi-mode at k=0. The thick solid lines show the dis- 
persions along [10] direction. The dispersions of the P-wave 
mode are compared very well with the EELS along [10] and 
[11] dispersions. See Ref. 3 for the detailed comparison. 



FIG. 5. Calculated intensities I (arb. unit) of the exciton 
transitions as function of k, normalized to Jo- The parameters 
are as same as those in Fig. 4. The thick solid lines show the 
intensities along [10] direction. The Pa-mode is not shown 
because it is dipole inactive and has no contribution to the 
intenisty for all k. The theory is compared well with the 
EELS along [10] and [11] directions. The intensities for the 
S-mode at finite k along [10] were observed in the EELS, see 
Ref. 3. 
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